Frequency dependent third cumulant of current in diffusive conductors 
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We calculate the frequency dispersion of the third cumulant of current in diffusive-metal contacts. 
The cumulant exhibits a dispersion at the inverse time of diffusion across the contact, which is 
typically much smaller than the inverse RC time. This dispersion is much more pronounced in the 
case of strong electron-electron scattering than in the case of purely elastic scattering because of a 
different symmetry of the relevant second-order correlation functions. 
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I. 



INTRODUCTION 



Measurements of nonequilibrium noise provide valu- 
able information about the properties of a system, which 
cannot be extracted from measurements of average quan- 
tities. For example, measurements of shot noise give the 
magnitude of the quasiparticle charge in the case of a 
tunnel contact and the effective temperature of electrons 
in the case of a diffusive contact J> Recently, Reulet et 
al^ performed first measurements of the third cumulant 
of current, which may give even more interesting infor- 
mation. For example, this cumulant is very sensitive to 
the presence of electron-electron scattering in a diffusive 
contact. Electron-electron scattering changes the shot 
noise in a diffusive contact only by several percent 
but it changes the third cumulant of current almost by 
an order of magnitude. 5 

Of special interest is the frequency dependence of the 
third cumulant. Very recently, it was shown that third 
cumulants of current in a chaotic cavity whose contacts 
have different transparencies may exhibit a frequency dis- 
persion much more complicated than that of the shot 
noise. Unlike the conventional shot noise that has a dis- 
persion only at the inverse RC time of the cavity^ the 
third cumulant of noise may also exhibit a dispersion at 
the inverse dwell time of an electron on the cavityi In 
most cases, this time is much longer than the RC time 
that describes charge relaxation in the cavity, and there- 
fore the corresponding dispersion takes place at exper- 
imentally accessible frequencies. This dispersion is due 
to slow fluctuations of the distribution function that do 
not violate electroneutrality and are akin to fluctuations 
of local temperature. These fluctuations do not directly 
contribute to the current and therefore are not seen in 
conventional noise, but they modulate the intensity of 
noise sources and therefore manifest themselves in higher 
correlations of current. 

Another important example of a system with a long 
dwell time is a diffusive contact. In this work we inves- 
tigate the frequency dependence of the third cumulant 
of a metallic diffusive wire. Like a chaotic cavity, it also 
has a long dwell time. In addition, the metallic diffusive 
wire is of interest because its impedance can be easily 
matched to that required by current experimental detec- 
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FIG. 1: The geometry of the system considered. White areas 
show the diffusive metal and dark areas, the insulator. 



tion schemes^ Furthermore, the measuring frequencies in 
this case are in the range where the frequency dispersion 
takes place for the system at hand. 

The zero-frequency third cumulant of current for 
a diffusive wire was first calculated by Lee, Levitov, 
and Yakovets& for non-interacting electrons in the zero- 
temperature limit. Recently this calculation has been 
extended to finite temperatures and to the case of strong 
electron-electron scattering by Gutman and Gefen£ In 
this paper we calculate the frequency dependence of this 
quantity both for the case of non-interacting electrons 
and for the hot-electron regime and show that the latter 
case is more convenient for the experimental observation 
of this effect, since the dispersion of the third cumulant 
is much stronger. 

We present a calculation based on the cascaded 
Boltzmann-Langevin approach^ In the Appendix, we 
also derive the full generating functional for the frequency 
dependent current fluctuations of a metallic wire both for 
the elastic and hot-electron regime based on the stochas- 
tic path- integral approach to full counting statistics. 10 
The third cumulant of current may be expressed in terms 
of functional derivatives of this functional. 
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II. MODEL AND BASIC EQUATIONS 

Consider a quasi-one-dimensional diffusive wire of 
length L and conductivity a. To explicitly describe its 
electric environment, the wire is chosen in the shape of 
a cylinder with a diameter 2rg and is embedded in a 
perfectly grounded medium, which is separated from the 
wire by a thin insulating film of thickness So and with 
a dielectric constant ea (see Fig. 1). All dimensions 
are assumed to be much larger than the elastic mean 
free path and the screening length in the metal. The 
electrodes are assumed to be perfect conductors, so the 
third cumulant of current is not affected by the external 
circuit We also restrict ourselves to sufficiently 

high voltages or temperatures, hence the quantum dis- 
persion of this quantity does not show up in the frequency 
range of interest^ We emphasize that despite the par- 
ticular choice of geometry, our results are valid for any 
quasi-one-dimensional diffusive contact. 

With the above assumptions, the noise of current 
may be described using the semiclassical Boltzmann- 
Langevin approach^- 5 The frequency dependence of shot 
noise in diffusive contacts with acco unt tak en of electri- 
cal screening was calculated in Refs. ITTITtL To calculate 
the frequency dependence of the third cumulant of cur- 
rent, we use the cascade extension of this approach^ The 
key idea of this extension is a large separation between 
the time scales describing the individual scattering events 
and the evolution of the distribution function of electrons 
in the contact. The resulting expressions may be also ob- 
tained by considering the corresponding stochastic path 
integraUS^ for the diffusive Boltzmann-Langevin equa- 
tion (see Appendix) . The cascade expansion corresponds 
to a systematic expansion of the saddle-point equations 
of this path integral in powers of the counting field. 

The quantity we are going to calculate is the Fourier 
transform of the third order current correlation function 
defined as 



ft(u>i,wa) = / d(h - h) J d(t 2 - h) 

xexp[ito 1 (t 1 -t 3 )+it02(t2-t 3 )](6I(t 1 )6I{t2)6I(t 3 )}.(l) 

The starting point for our calculations is the stochastic 
diffusion equation for the fluctuations Sf(e, r) of the dis- 
tribution function f(e,r) 







dt 



Of 



tt-DV 2 ) 5f-6I ee = -eS^-VSF lm P-8F ee , (2) 



ds 



where D is the diffusion coefficient, SI ee is the linearized 
electron-electron collision integral, and <5F lmp and 5F ee 
are random extraneous sources associated with electron- 
impurity and electron-electron scattering. This equa- 
tion is obtained from the standard Boltzmann-Langevin 
equation by defining the electron energy as e = p 2 /2m + 
e(f>(r , t) — ep and isolating the isotropic part of the distri- 
bution function in the momentum space. The fluctuation 
of the electric potential 5<p that appears in this equation 



should be calculated self-consistently from the Poisson 
equation 

V 2 S(j) = -inSp, (3) 
where the fluctuation of charge density Sp is given by 

Sp = eNp (^J deSf{e) + e5(^j (4) 

and where Np is the Fermi density of states. In the case 
of a quasi-one-dimensional contact, a solution of Eqs. (J2J 
- (HJ is of the formil 

5<j)(x,uj) = J— (V 2 +iLuRC/L 2 Y 1 (5) 
boa 



xiz I d 2 r ± 6f x xt (v), 



d_ 

dx _ 

where x is the coordinate along the contact, a — e 2 NpD 
is the conductivity of the metal, So = tttq is the cross sec- 
tion area of the contact, C — Ledro/25o and R = L/irr^cr 
are the capacitance and the resistance of the contact, and 



eN F / de5F mip . 



(6) 



A fluctuation of the total current density is given by 

Sj = 5r xt - tVJ0, (7) 

and a fluctuation of the total current at the left end of 
the contact thus equals 



81 = a 



dS(f>(x, lj) 



dx 



-L/2 



(8) 



Making use of the correlation function of extraneous 
sources 



D 



(5Fr P (e,r)<5Fr p ( £ ',r')) w - 2— S(r - r')5(e - e' 



x6 a pf (e,r)[l 



N F 



(9) 



one easily obtains the second-order correlation function 
for the fluctuations of the current as a functional of the 
distribution function 

Consider now the third cumulant of current. As the 
direct contribution to this quantity from the third cumu- 
lant of extraneous sources is negligibly small in a diffusive 
metalfS this quantity is dominated by an indirect contri- 
bution of the second cumulant of these sources, which 
results from the modulation of their intensity by fluctu- 
ations of the distribution function. It may be written in 
the form 



(SI(h)6I(t 2 )SI(t 3 )) = Pi 23 {A 123 }, 



(10) 



where 



A 



123 



dt / de 



Sf(e,r,t) 
x{S.f(e,r,t)SI(t 3 )) 



(11) 
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and V123 denotes a summation over all inequivalent per- 
mutations of indices (123). 

Equations © and (JSJ suggest that the second cumu- 
lant of current exhibits a dispersion at frequencies of the 
order of (i?C) _1 . Typically such high frequencies are be- 
yond the experimentally accessible range. Therefore in 
what follows we will assume that the frequencies Wi, u>2, 
and L03 are much smaller than (i?C)~ 1 . Hence the pile- 
up of charge in the contact may be neglected and the 
fluctuation of current may be considered as coordinate 
independent. In this case, 

{6I{uji)5I{lu 2 )) = 4vr(5(wi + ^{RLy 1 (12) 
x I dx j def(e,x)[l - f(e,x)] 



and the only possible dispersion in Eq. (|11|) is due to the 
dynamics of a fluctuation 8f, so that the expression for 
the third cumulant assumes the form 

P 3 (u>i, w 2 ) = P(wi) + P(wa) + P{~ui - ^2), 

L/2 



P{uj) = — J dxlds[l-2f(e,x)}{6f(e,x)5I} u .(13) 

-L/2 

The quantity P(tu) has to be calculated in different ways 
for the case of purely elastic scattering and for the hot- 
electron regime. 
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FIG. 2: The coordinate dependence of the zero-frequency 
correlator (5<j)(x)SI) ul normalized to eV for the elastic case 
(solid line) and the hot-electron regime (dashed line). 



and 



(SU(s, x)SI) u = iuj (V 2 - iu/D) 1 
~df(e,x) 



ds 



-e{8(j>(x)8I} u 



(17) 



At low frequencies, one easily obtains from Eqs. 10 and 
JHl that 



III. 



PURELY ELASTIC SCATTERING 



For purely elastic scattering, SI ee and SF ee in Eq. @ 
vanish, and a fluctuation of the distribution function Sf 
may be presented as a sum of a part induced directly by 
an extraneous source 



6f F (e,x,u) = (V 2 + ilo/D) 1 



dx 



(14) 



and a part induced by fluctuations of the electric poten- 
tial 



£, x, u>) — ~ iu> (V 2 + iu>/D) 



df(s,x) 



{x,u) 



(15) 

The existence of the term l|15f) indicates that the dy- 
namics of charged electrons differs from the dynamics of 
neutral particles even at frequencies much smaller than 
(RC)-\ 

By multiplying these equations with the fluctuation of 
current (|SJ) and making use of the correlation function 
@, we obtain 



(Sf F (e,x)SI) u 



-2|(V 2 



} /D)- 1 



x—{f(e,x)[l-f(e,x)]} (16) 



{84>(x)8I) u - I (V 2 )" 1 



lxI dS ^~ f \ 



(18) 



Using the well-known expression for the average distri- 
bution function 



f(s,x) = [- + -)f (e+eV/2)+ 



f (e-eV/2), 
(19) 

where fo is the equilibrium Fermi distribution and V is 
the voltage drop across the contact, we obtain 



(§4>(x)si) u 



1 X 

6L 



1 



eV coth 



eV 
2T 



2T 



(20) 

The correlator (5cj)(x)5I) vanishes at V = and is an odd 
function of x at nonzero V(see Fig.|2J). Upon inverting 
the operator (V 2 + iuj/D) in Eqs. JTHJ) and l(T7|l and 
performing the spatial integration in Eq. Ijl3[) . we arrive 
at an expression for P{ui) in terms of q u — (zw/D) 1 / 2 , 
which is our final goal. Because of its length, we give 
here only its low-temperature and low-voltage limits 
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4e 2 U, 



qMqtL' + 30) sinh(^L) 



3 R 

8(qlL 2 + 6) cosh^L) + 2qlL 2 
/[&L B swh(q u L)] 



•48] 



(21) 



4 




0.06 



0.05 



0.04 



40 , - _ 60 

cox D 



(M 

^) 

n 0.03 
E 

0.02 
0.01 





— elastic, T=0 




- - inelastic, T=0 




elastic, eV->0 




— inelastic, eV->0 ' 


• : ' * '■. 




;V \ \ 




i V ' 




■ ' \ ■ 




■i ' \ \ \ 




V N V\ 

i 




!< s % •>•••. 








w 








? 




1 







20 



40 cox D 60 



100 



FIG. 3: The real part of the ratio P(to)/ (e 2 I) versus normal- 
ized frequency lotd (td = L 2 / D) for purely elastic scatter- 
ing at eV 2> T (solid line), hot-electron regime at eV 3> T 
(dashed line), purely elastic scattering at eV <C T (dotted 
line), and hot-electron regime at eV <C T (dash-dotted line). 
Inset: charged and uncharged fluctuations of the distribution 
function /. Charged fluctuations have a short relaxation time 
and contribute to fluctuations of current SI. Uncharged fluc- 
tuations do not contribute to SI directly but affect the inten- 
sity of noise sources. They decay only via slow diffusion and 
result in the low- frequency dispersion of the third cumulant. 



and 



4 e 2 V 2 cosh(q w £) - q u L sinh(g^L) — 2 



3 R 



q 3 L 3 sinh(q u L) 



(22) 



At oj — these expressions give — (l/45)e 2 V/i? 
and — (l/9)e 2 V/ R, which corresponds to P3(0,0) = 
-{l/lh)e 2 V/R and P 3 (0,0) = -(l/3)e 2 V/R. These 
zero-frequency results are in agreement with Ref. 00 . 
At finite frequency, Eqs. I|21l) and i|22|) become complex- 
valued and tend to zero as i/ui at lo — » oo. The real and 
imaginary parts of P(lu) are shown in Figs. and 01 



IV. 



HOT-ELECTRON LIMIT 



Consider now the limit of strong electron-electron in- 
teraction. In this case, the distribution function may 
be assumed to have a Fermi shape with a coordinate- 
dependent temperature T e (x) and electric potential <p(x) 



1 + cxp 



T e {x) 



(23) 



If the frequency to is smaller than the rate of electron- 
electron collisions, a fluctuation Sf can be expressed in 
terms of fluctuations of these quantities 



5f(e,r,w) 



df(s,r) 



df(e,r) 
dT e 



6T e 



(24) 



FIG. 4: The imaginary part of the ratio P(u>)/(e 2 I) ver- 
sus normalized frequency lotd (td = L 2 /D ) for purely elas- 
tic scattering at eV 3> T (solid line), hot-electron regime 
at eV S> T (dashed line), and purely elastic scattering at 
eV <C T (dotted line), and hot-electron regime at eV <C T 
(dash-dotted line). 



A substitution of Eq. I|24|) into Eq. I|13|) and integration 
over the energy readily gives 



L/2 



(25) 



-L/2 



To calculate the correlator in Eq. (|25|) , we have to obtain 
a Langevin-type equation for 5T e . To this end, we multi- 
ply Eq. J2| by e and integrate it over e, like it was done 
when deriving the equation of heat balance in Refs. 1.31 ll 
This gives 



d 



—T e 6T,, I - DY 



2t e 2 



deeYSF 1 



(26) 



Multiplying Eqs. H2bj) and |J8J| and averaging the product 
with the help of Eq. © results in an equation for the 
correlation function (5T e 5I)u of the form 



*"-<b) 



-T e {5T e {x)8I) u 



-V 2 [e^S^SI)^ + J deef(l - /). (27) 

The integral over the energy on the right-hand side of 
Eq. (|27|l equals e(f)T e , and making use of Eq. (|I8(I . one 
easily obtains the solution of Eq. 1)2 7|1 in a symbolic form 



6 e 2 1 



(ST e (x)SI) u = - --(V^ 

7T Z L lp 



f d(<f>T e 
\ dx 



0(V 2 



iuj/D)' 

dTe 
dx 



(28) 
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where the operators V 2 and V 2 — iui / D are inverted with 
zero boundary conditions. According to Ref. 0, the mean 
potential is given by <f>{x) = Vx/L and the mean effective 
temperature by 



T e (x) = 



HP* 



3 

— 1 





x 2 Y 


(i- 


L 2 ). 



1/2 



(29) 



As <p and T e are odd and even functions of coordinate, 
the resulting correlator is an even function of x. In the 
zero-frequency limit it is given by 



infinite series 



Phot{u) = — — 

k=0 



7T R 



Jo( 7 rfc + 7 r/2)[J 1 (7rfc + 7 r/2)-(-l) fc ] 

(2k + l)[TT 2 (2k + l) 2 + iluL 2 /D] ' { } 

where Jo and J\ are Bessel functions of order and 1. 
V. DISCUSSION 



(5T e (x)SI) u 



1 



^a 2 -4x 2 /L 2 
+6 



(a 2 _ 1} 3/2 _ 3 * arcsin ( |£ \ 



a arcsm 



a; 
L 2 



-,(30) 



where a 2 = 1 + (47r 2 /3)T 2 /(eF) 2 . This correlation func- 
tion vanishes at V = and is negative at positive volt- 
ages, i.e. if the increment of total current through the 
contact is negative. This fact has a very simple physical 
meaning: an increase in the total current results in an 
increase of the Joule heating and hence an increase of 
the temperature. 

An integration of Eq. I|28|l with respect to x gives the 
low-frequency third cumulant for arbitrary temperatures 
and voltages in a form 



P 3 (0,0) = 4^ 
n z R 



^-H a2+ i (5a 2 -2)v^T 



4 -2 

-a arcsm 



(31) 



Its limiting values P 3 (0,0) = -{3/ir 2 )e 2 V/ R at eV < T 
andP 3 (0,0) = -(8/ir 2 -9/16)e 2 V/ R at eV > T coincide 
with the results of Gutman and GefenA 

In the case of nonzero frequencies, it is possible to 
obtain analytical results only in the limiting cases of 
eV <C T and eV S> T. In the high-temperature limit, 
one may set T e (x) = T in Eq. (J2BJ, hence <j>(x) is 
the only coordinate-dependent quantity, and the term in 
curly brackets equals VT/L. Then the diffusion operator 
is easily inverted and the integration over x gives 



Phot{u) = - 



12 e 2 V 1 



1 tanh 



2 



(32) 

In the zero-temperature limit, the term in curly brack- 
ets in Eq. (|28|l may be expanded in a Fourier series in 
cos[(2n + \)-kx/L] and the operator V 2 — iu/D is eas- 
ily inverted. The final result is obtained as a sum of an 



Our results for the non-interacting regime (Eq. 12 II and 
Eq. I22fl and for the hot-electron regime (Eq. [221 and Eq. 
|33Jl are displayed in the figures and ^ It is clearly seen 
that both real and imaginary parts of the third cumu- 
lant have the most pronounced dispersion in the case of 
a high temperature or for strong electron-electron scat- 
tering, i.e. when the local distribution function has a 
nearly Fermian shape. This unexpected result is in a 
sharp contrast with the dispersion of quantum noised 
which results from sharp singularities in the energy de- 
pendence of the distribution function. 

Mathematically, the different shape of the frequency 
dependence for purely elastic scattering and interact- 
ing electrons at high voltages can be explained as fol- 
lows. In the case of hot electrons, both S(5I 2 ) /ST e (x) 
and {5T e {x)8I) UJ are even functions of the coordinate 
x (measured from the middle of the contact). On the 
contrary, for purely elastic scattering at T = both 
8(8I 2 ) /5f(e,x) and {5f(e,x)8I) UJ are odd functions of x. 
The functions acted upon by the inverse diffusion oper- 
ator (V 2 — iuj/D)^ 1 in Eqs. (|2*S|) and ljTH are also even 
and odd, respectively. At low frequencies, the inverse 
diffusion operator is essentially nonlocal in space and ap- 
plying it results in an effective averaging of the argument 
on the scale of the order of L. In the elastic case, this 
averaging involves both negative and positive values, and 
this is why the elastic third cumulant is suppressed at low 
frequencies as compared to the hot-electron value. How- 
ever at high frequencies, the inverse diffusion operator 
becomes almost local in space, therefore there is no aver- 
aging of negative and positive values and both cumulants 
become of nearly the same magnitude. This absence of 
spatial averaging partially compensates for the increas- 
ing frequency and makes the frequency dependence of 
the "elastic" cumulant more flat. Therefore the differ- 
ent shape of the frequency dependence in the elastic and 
hot-electron limits may be traced back to the different 
symmetry of relevant second-order correlation functions. 



VI. 



SUMMARY 



In summary, we have shown that diffusive contacts ex- 
hibit a nontrivial internal dynamics even at frequencies 
much smaller than the inverse charge-relaxation time. 
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Though this dynamics is not affected by the electric en- 
vironment of the contact, it differs from the dynamics 
of charge-neutral particles and manifests itself as a low- 
frequency dispersion of the third cumulant of current. 
In view of the fact that both dynamic conductance and 
shot noise of metallic conductors depend only on the RC 
time, this frequency dispersion of the third cumulant on 
the scale of the dwell time is a very interesting result. 
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APPENDIX A: STOCHASTIC PATH INTEGRAL 
REPRESENTATION 

In this appendix, we derive a stochastic path integral 
representationi2ii& for the full current statistics of diffu- 
sive conductors. This representation then serves to cal- 
culate the dispersion of the third order current correla- 
tion function. We give first a detailed derivation of the 
case without electron-electron scattering. For simplicity, 
we restrict ourselves to quasi-one-dimensional wires and 
consider only slow dynamics due to the longitudinal dif- 
fusion modes. In this case, the kinetic transport equation 
for the fluctuating contribution Sf(e, x, t) to the electron 
occupation function / = f(e,x) + 5f(e,x,t) is given by 
the one-dimensional equivalent of Eq. J5J). The mean 
distribution function / is defined in Eq. I|19|) . Since we 
restricted ourselves to the regime of slow diffusive dy- 
namics, we may assume charge neutrality d{8j)/dx = 
where the fluctuations of the electrical current density 
are given by Eqs. © and J7J. Finally, we make use of 
the facl£ that the extraneous sources of noise are Gaus- 
sian random variables delta-correlated in time, space and 
energy and described by Eq. |(SJ|. 

In order to construct a path integral representation of 
this Boltzmann-Langevin equation, we define the proba- 
bility functional P which gives the probability to find a 
certain realization of extraneous currents SF lmp (£,x,t). 
This functional may be written as a path integral 

P [SF" n P] = 
J Vr)cxp{f dtdxde{~ir,8F lmp + H(r]))} (Al) 

taking Fourier transforms of the generating function 
D 



H 



N F 



f(i-fW 



(A2) 



Since the Fourier transforms are independently taken for 
each point in space, time and energy, Eq. (|A1() indeed 
characterizes white noise. The electron occupation func- 
tion / = f+5f is considered to be a slowly changing vari- 
able which modulates the instantaneous noise intensity. 
Its evolution is determined by the kinetic equation (J2J) 
which we represent as a delta functional expressed by a 
path integral 



Sf - DSf" + e^Scj) + (SF" np Y 



= / VXexpi 



de 



dtdxdeX [Sf 



(A3) 



DSf" 



de 



(SF 



imp\f 1 



where the prime stands for d/dx. The dynamics of the 
potential fluctuations S<p can be expressed by a second 
delta functional which enforces charge neutrality 



S [c 



zN F J de(SF imp )'} =JV£, (A4) 



exp {i J dtdx£, (a5(f>" ~ eN F J de(5F lmp )')} . 

The fields A and £ can be understood as Lagrange mul- 
tipliers. Combining Eqs. (|A1|) . (|A3|> and l)A4|l . we con- 
struct the probability Pt to find a certain realization of 
extraneous currents SF" np under the constraint of cur- 
rent conservation and charge neutrality 

P t [5F imp ] = J VS<pVSfS[, . . ]S[. . . ]P[SF imp ]. (A5) 

We are now in a position to calculate the generat- 
ing functional S[i\] of current fluctuations SI — Sj(x = 
—L/2) at the left contact x = —L/2 

e Sc,[lxl = J VSF mp P t [SF mp ]exp^i J dt X SI^ . (A6) 

This equation may be considerably simplified. In a first 
step, we can integrate out the the extraneous currents 
SF lmp as well as the field r\ introduced in Eq. IIAlll . 
We are then left with four functional integrations over 
Sf, 5<p, A and £. In a second step, we evaluate this inte- 
grations in the saddle point approximation. As the diffu- 
sive conductor is essentially semiclassical, the corrections 
to the saddle point action are smalUS. After rescaling 
A i— » eNpX, we are left with the generating functional 



S e i\x, A, t Sf, 5^] = f dtdxfasr + f 



de 



fA<5f+ 

'TJE^)]} ( A7 ) 
which has to be evaluated at the saddle point given by 



Tf(l- /)(A' + 2 -eN F X[Sf + eM. 



SS e l 

sx 
ss el 

S(Sf) 



= 0, 



0, 



SS e i 
St 

ss el 

5(S<t>) 



= 0, 
= 0. 



(A8) 
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Note that we performed a complex continuation iX i— » 
A, i£ i— > £ and ix i— » %. We are therefore left with 
purely real quantities. The saddle point equations are 
supplemented with boundary conditions: the three fields 
Sf, 6cj>, A vanish at both boundaries. The external count- 
ing field x is incorporated into the boundary conditions 
for £ 



£(-L/2) = X , £(V2)=0. 



(A9) 



The frequency dependent third cumulant under consid- 
eration in this paper is obtained from the third functional 
derivative 

S 3 S el 



{6I(ui)5I(u 2 )SI(lu 3 )) 



<5x(wi)£x(^2)<5x(w3) 



(AlO) 



x=o 



We calculate this cumulant by a systematic expansion of 
action JXJJ, saddle point equations <| A8|> and fields 



Sf = 5f x + 5f 2 + 



= 5(f)! + 



(All) 



in orders of the external field x- It can be straightfor- 
wardly shown by inserting saddle point equations (|A8|) 
back into the action (|A7(1 that the third order contribu- 
tion to the action has the form 



SeiAx] =o I dt J dx I de(l-2f)5h(^Y. (A12) 

It remains to solve the saddle point equations to first 
order in x- For the Lagrange multipliers we find 

K +f" = 0, thus Ai = 0,6 = ~x(x/L - 1/2). (A13) 

Their dynamics is trivial, since they follow instanta- 
neously the external field x- The interesting dispersion 
effect stems from the saddle point equation for the oc- 
cupation function which takes the form of an inhomoge- 
neous diffusion equation 



Sfi - DSfi' 



(A14) 



~2eD (/ (1 - fo)^)' ~ 2eD^ J de (/ (1 - /„)£)' . 

This diffusion equation has a very appealing interpreta- 
tion: When we integrate the equation over energy, we find 
that the two source terms on the right hand side cancel. 
The left side becomes a homogeneous diffusion equation 
for the fluctuations of the charge density 8p\ = e J deSfi 
which has the trivial solution Spi = 0. This is nothing 
than the charge neutrality which we demanded in the 
beginning of this section. The second source term which 
is due to variations of the electrostatic potential thus 
compensates the first term in such a way that all fluctu- 
ations of the occupation function Sfi are charge neutral. 
We decompose the total variation Sfi = Sf[ + Sff into 
a contribution Sfi due to free fluctuations of the occu- 
pation function and a contribution Sff due to potential 
fluctuations. Using the identity 5<j>i — DS<fi", we solve 
the diffusion equation by inverting the diffusion operator 
and find 



6f[(cj) = (5f F (e,x)SI) u 
6tf(w) = {6U(e t x)6I) u 



D ' 

xH 

D ' 



(A15) 
(A16) 



The two correlators are defined in Eq. (|16fl and Eq. i|17|) 
respectively. They can now be inserted into Eq. I|A12(I to 
obtain the third-order correlation function Ijl3|l . We thus 
derived the cascade rules applied in the main body of this 
paper from the stochastic path-integral formalism. 

The derivation of an action describing the hot electron 
regime requires only a minor additional effort. Here we 
cite directly the result which has been derived by one 
of the authors for the zero frequency limit in a different 



context 



Shot = / dt f dxl -N F XT e f e 



+a(C X')A[ ) -a(£ X')B ( ) } (A17) 



T 



In this action, we introduced the local electron temper- 
ature T e (x,t) — T e (x) + STe(x,t) and the local electro- 
static potential 4>(x, t) — <p{x) + 5<f>(x, t). The boundary 
conditions for these fields are the potentials and temper- 
atures of the left and right reservoirs. As for the case 
of non-interacting electrons, the Lagrange multiplier £ 
ensures charge neutrality and obeys the boundary con- 
dition (|A9|) . The field A is linked to the conservation of 
energy current and is zero at the boundaries. The ma- 
trix A describes the local noise created by the extraneous 
sources of noise SF lmp 



A = T R 



1 



(0 2 +^ 2 T e 2 /3e 2 ) 



(A18) 



The second matrix B is the linear response tensor 



B = 





7r 2 T e /3e 2 



(A19) 



In complete analogy to the derivation of Eq. (|A1 21) , we 
may again collect all third order terms which contribute 
to the action l|A17(l and find 

S h ot,z[x] =o J dt f dx5T ejl (Z[) 2 . (A20) 



where the variation ST e i can be identified with 



ST e ,i = (ST e (x)SI) u 



(A21) 



(see Eq. (|28H 'I. The total third order contribution 
Shot,3[x] therefore corresponds exactly to Eq. 

The main results of this appendix are the dynamic gen- 
erating functionals Eq. (|A7|) and Eq. (|A17|) . These func- 
tional permit (in principle) the calculation of cumulants 
of arbitrary order. 
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